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ABSTRACT 

We report seven successful observations of the astrometric binary GJ 164 AB system with aperture 
masking interferometry. The companion, with a near infrared contrast of 5:1 was detected beyond 
the formal diffraction limit. Combined with astrometric observations from the literature, these ob- 
servations fix the parallax of the system, and allow a model-independent mass determination of both 
components. We find the mass of GJ 164B to be 0.086 ± O.OOTMq. An infrared spectroscopic study 
of a sample of M-Dwarfs outlines a method for calibrating metallicity of M-Dwarfs. Results from the 
newly commissionned TripleSpec spectrograph reveal that the GJ 164 system is at least of Solar metal- 
licity. Models are not consistent with color and mass, requiring a very young age to accommodate a 
secondary too luminous, a scenario ruled out by the kinematics. 

Subject headings: binary, stars: luminosity functions, mass functions, techniques: AO, interferometric 
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1. INTRODUCTION 

The development of a self-consistent theory of the in- 
ternal structure of stars and the construction of models 
which trace evolutionary behaviour are major achieve- 
ments of modern astrophysics. Both theory and models 
have however long been centered on intermediate and 
high mass stars, with little consideration of objects with 
masses below 0.6 Mq. The realization that the Solar 
neigh borhood is overwhelmingly dominated by low mass 
stars (|van de KainpllQTltlHenrvjflQQS.). a nd the discovery 
of brown dwarfs ( Nakaiima et al.lll995f) have brought a 
lot more attention to the lower end of the main sequence 
and what lies below. 

The complex physical nature of these very low mass 
objects makes their modeling intricate. Boundaries be- 
tween stars, brown dwarfs and planets, as well as the- 
oretical relations that predict the luminosity and tem- 
perature as functions of age, mass and metallicity are 
still largely untested in the relevant range of parameters. 
Correct calibration of these models is of extreme impor- 
tance, since they are now extrapolated to estimate the 
mass of b rown dwarfs (iBurrows et al.ll2006 h and giant ex- 
oplanets (|Baraffe et al.l 120031 ) from their luminosity and 
age. 

The observations required to challenge and improve 
the models are dynamical mass measurements of mul- 
tiple star systems, combined with accurate photome- 
try and dista nce determinatio n. Besides gravitational 
microlensing (jPaczvnskil Il986f ). the observation of bi- 
nary systems and the use of Newtonian orbital dynam- 
ics provide the only method of directly measuring ac- 
curate stellar masses. With maximum sensitivity at 
separations greater than 1 arcsec, conventional Adap- 
tive Optics (AO) observations at Palomar are limited 
to the most nearby objects or the consequent orbital 
periods become too long to lead to dynamical masses. 
We have therefore been using aperture masking inter- 
ferometry ( Tuthill et al.l l 2000D in conjunction with AO 
(|Llovd et al.ii2006ll . The precision calibration of the data 
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achieved with this observing mode indeed leads to reli- 
able results up to and beyond the formal diffraction limit 
(super resolution), and ve ry precise photornetry. Com- 
bined with rad ial velocity dMartinache et al.l [ 2007| ) or as- 
tromctry (Irel and et al.l 120081 ). aperture masking inter- 
ferometry has so far provided some of the most precise 
(dynamical) masses of objects below 0.1 Mq. 

2. CHARACTERIZATION OF THE ORBIT 

We report here aperture masking interferometry ob- 
servations of the astrometric binary GJ 164 AB. GJ 164 
(aka LHS 1642, G 175-19 or Ross 28) is a high proper 
motion star, catalogued as a M4.5 Dwarf (jReid et all 
1 1 995( 1 ■ 11-13 pc distant. From the 2MASS catalog, its 
apparent near infrared magnitudes are J = 8.773 ±0.032, 
H = 8.248 ± 0.030, and K = 7.915 ib 0.016. The discov- 
ery of GJ 164 B was reported by iPravdo etHII (12004') . 
as a p art of the STEPS program of iPravdo fc ShaklanI 
(|1996l ). The combination of the astrometry data from 
the discovery paper, combined with the aperture mask- 
ing interferometry data reported in this paper provides 
a complete dynamical characterization of the system, re- 
sulting in masses that are independent of the use of a 
model or a mass-luminosity relation. GJ 164 is an espe- 
cially interesting target for the mass of the B component 
lies very close to the substellar limit. 

2.1. Aperture Masking Interferometry 

Measuring the dynamical masses of a low-mass binary 
system requires some patience, for the observations have 
to cover at the very least a significant fraction of the 
orbital period. Because of its low mass (< 0.6 Mq), 
a relatively short period (< 10 yr) binary will have a 
semimajor axis smaller than 4 AU. Even for objects as 
close as 10 parsecs, the angular separation won't exceed 
0.4 arcsecs, a performance theoretically within the grasp 
of a 5-10 meter class telescope equipped with AO. Yet 
in practice, the residual quasi-static speckles in the AO 
PSF halo seriously limit the capability of detecting a 
companion at angular separation smaller than 2-4 \/D 
([Racine et al.l 11999) . Aperture masking interferometry 
with AO adresses this issue. 

A non redundant mask discards most of the pupil by 
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Fig. 1. — Calibrated closure phases observed on JD 2453779.71 
with a Kg filter. A binary star model (light gray) with a 86 mas 
separation, a position angle of 102° and a contrast of 5:1 satisac- 
torily matches the observations (dark gray). 

sampling a few spatia l frequencies only (Readhea d et al.l 
[198 8: Nakaii ma et al.l[r98 9). One admittedly loses most 
of the light (the 9-hole mask used for this works trans- 
mits approximatively 15 %) but rejects all atmospheric 
noise, as well as internal aberrations and non common 
path errors, which i n turn, dra matically increases the 
signal-to-noise ratio (jTuthill et al. 20001. Used in con- 
junction, aperture masking interferometry and AO pro- 
vide stable fringes, and enable long integration times, 
therefore making faint target accessible. The only limit 
is whether the AO system can lock on the system. Thus, 
any target that can be observed with AO can also be 
observed with aperture maskin g interfero metry. With 
a visible magnitude V = 13.5 (|Weislll996h . GJ 164 can 
fairly easily be observed with aperture masking. If cal- 
ibration is sufficiently accurate, the interferometry en- 
ables super-resolution, ie. the detection of structures on 
a target at a scale smaller than X/D. The data presented 
in Tables [Hand [2] shows that the companion to GJ 164 A 
is detec ted on multiple occa sions in that super-resolution 
regime. iTuthill et al.l (|2006D provide a recent description 
of the general principles and performance s of a perture 
masking interferometry while iLlovd et al.l (|2006D detail 
the experiment undertaken at Palomar. 

The AO observa tions performed with PHARO 
(jHayward et al.ll2001f ) at Palomar as well as with NIRC2 
at Keck II span the range December 2003 to December 
2006. The companion was successfully detected all seven 
runs. The data are calibrated with observations of one or 
several stars of similar brightness for both the AO wave 
front sensor and the science camera to ensure compara- 
ble wave front correction and signal-to-noise ratio. The 
final dataproduct of our custom software pipeline writ- 
ten in IDL, is a coll ection of calibrated closure phases 
([Baldwin et al.lll986D . This very robust observable re- 
jects both atmospheric noise and calibration errors of the 
wave front sensor. The 9 hole mask provides (|)=84 pos- 
sible closure triangles, (f) — 28 of which are theoretically 
independent. 

The instantaneous geometry of a binary star is mod- 
eled by a set of three parameters: contrast ratio c, angu- 
lar separation a and position angle (PA) 9. The example 
of one such model fitted to a set of closure phases is given 
in Figure [TJ The aperture masking data processing soft- 
ware typically uses several hundreds of frames which are 



TABLE 1 

APERTURE MASKING MEASUREMENTS AT 
PALOMAR AND KECK: angular separation and 
POSITION angle of GJ 164 B. 



Julian Date 


Band 


Telescope 


Sep. 


PA 


(-2,450,000) 










3004.85 


H 


Palomar 


82.8 ± 1.6 


111.1 ± 0.9 


3632.96 


H 


Palomar 


53.3 ± 1.3 


149.1 ± 1.2 


3723.66 


H 


Palomar 


80.1 ± 3.3 


116.1 ± 2.1 


3779.71 


H, Ks 


Palomar 


86.1 ± 2.0 


102 1 -1-0 9 


3957 62 


H 


Keck 


49.3 ± 0.5 


47.0 ± 0.6 


401 8 98 

rtVJ J.O. i/O 


Ks 


Palomar 


39.4 ± 1.3 


■^44 8 -1- 2 6 


4078.79 


Ks 


Palomar 


57.3 ± 2.2 


294.1 ± 2.8 






TABLE 2 








GJ 164 PHOTOMETRY 






Julian Date Filter 


Contrast 






(-2,450,000) 








3004.85 


H 


5.97 ± 0.32 






3632.96 


H 


5.77 ± 0.46 






3723.67 


H 


5.78 ± 0.63 






3779.71 


H 


5.53 ± 0.60 








Ks 


4.89 ± 0.19 






3957.62 


H 


5.42 ± 0.03 






4018.98 


Ks 


8.66 ± 4.81 






4078.79 


Ks 


5.67 ± 1.17 





averaged to produce the closure phases. Due to the cen- 
tral limit theorem, the distribution of closure phases {0^} 
can be considered Gaussian, with uncertainties oi given 
by the standard deviation to the mean. The best fit mini- 
mizes the traditionally used goodness-of-fit parameter ^ 
defined as 

2 _ - MO',S,cjf 

i * 

where fi designates the model. A systematic error term 
is added to the closure phase dispersion Ui to achieve a 
reduced xl = 1. The likelihood of the parameters given 
the set of closure phases {cjij} is related to the by 

i(a,0,c|{0,})(xexp(-xV2). (2) 

Normalized, L is the full joint probability density func- 
tion for all three parameters. The confidence interval of 
one individual parameter is calculated by integrating out 
the two others. For example, the marginal probability 
density function for the angular separation is calculated 
by integrating out the position angle and the contrast 
ratio: 

p{a) ^ j dedcL{a,e,c). (3) 

The results of this analysis: relative position and pho- 
tometry are gathered in Tables [1] and [5| The more or 
less favorable seeing conditions explain the variable con- 
fidence interval sizes. We achieve a precision of a few 
milliarcseconds at separations as low as 44 mas in K5 
band (0.5 A/d) on good nights at Palomar. 

2.2. Astrometry 

G J 164 was discovered as a b i nary i n an astromet- 
ric survey by iPravdo fc Shaklar] ([199(3( ): iPravdo et all 
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TABLE 3 
GJ 164 ASTROMETRY 



ian Date 
450,000) 



AR.A. A Decl. 



D801.5 
D801.5 
1088. 3 
L189.0 
L189.0 
L189.0 
L436.5 
L494.0 
2657.5 
2889.5 
2889.5 
2889.5 



0.00 
0.00 
-147.91 
-328.82 
-324.82 
-325.32 
-457.34 
-563.37 
-1639.04 
-1745.18 
-1745.18 
-1745.88 

-9014 .'^1 



0.00 
-2.80 
-646.96 
-874.15 
-874.05 
-871.65 
-1409.82 
-1519.60 
-4132.08 
-4629.59 
-4627.29 
-4625.69 



100 



o 
n 




(-2. 



3030.5 -2014.61 -4957.45 
3051.5 -2044.87 -5016.85 



(|2004[ ). While the interferometry resolves the binary and 
provides its instant geometric configuration, the astrom- 
etry records the position of a star's photocenter, mea- 
sured relative to several more distant stars. Repeated 
observations reveal the star's proper motion as well as its 
parallax. Once these two effects substracted from the as- 
trometric signal, any residual wobble of the photocenter 
can be attributed to the presence of one or more com- 
panions, massive enough to cause an observable shift in 
position. The orbit of the companions can then be fully 
characterized. 

Astrometric data were extracted from the figures of 
iPravdo et al.l (|2004l ) then added the motion due to the 
proper motion and the parallax given in the paper's ta- 
bles. These extracted values are gathered in Table |3l We 
initially assume a 2.0 mas uncertainty for each individ- 
ual measurement, as suggested by t he re siduals of the 
analysis performed bv iPravdo et al.l (|2004f ). 

The astrometric residual wobble however remains a 
function of the fractional light of the secondary, /3 = 
L2/{Li + L2), with Li and L2 respectively denoting the 
luminosity of the primary and the secondary, as well as 
the fractional mass of the secondary / = 7712/ (mi -I- m2). 
In the limit at which the light of the secondary is negli- 
gible, the photocentric orbit is identical to the Keplerian 
orbit of the primary around the actual center of mass of 
the system. If the luminosity of the secondary is not so 
negligible, then the photocentric semimajor axis a is re- 
duced proportionally to (3. The ratio of the photocentric 
orbit a, and the Keplerian orbit a is (jMcCarthv et al.l 
fl988il : 

a/a^f~ 13. (4) 

A relatively small-mass, non-luminous secondary, and 
a relatively large-mass, luminous secondary will therefore 
have indentical astrometric signatures. The data pre- 
sented in Section 12.11 however rules out the latter possi- 
bility. Indeed, the aperture masking interferometry data 
gathered in Table [2] provides the H and Kg band contrast 
ratios, which are turned into the following differences of 
magnitude: 

AiJ= 1.835 ±0.006 
AisTs^^ 1.721 ±0.097. 
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Fig. 2. — GJ 164 B ORBIT. The measurements and associ- 
ated uncertainties of the Aperture Masking data are represented 
by dark gray rectangles. Light gray triangles represent the STEPS 
astrometry data. The relative apparent orbit (curve) is clockwise. 
Table |4] provides the corresponding orbital elements. The almost 
horizontal straight line represents the line of nodes. 

These contrast ratios can be used to decompose the 
2MASS combined magnitudes of the binary into ap- 
parent magnitudes for individual components: H = 
8.432 ± 0.030 and K = 8.117 ± 0.017 for GJ 164A, 
H = 10.267 ± 0.030 and if = 9.834 ± 0.038 for GJ 164B. 

The corresponding (H-K) color indices: {H — K)a = 
0.315±0.034and {H~K)b = 0.433±0.048 alone, provide 
a parallax-independent observable. While the {H — K) 
color index of the primary is well compatible with a M4.5 
spectral type, the secondary is v ery unlikely to be any 
earlier than M8 (see for instance iKoornneea (|1983D ). A 
well characterized nearby equivalent is the star VB 10 
from the 8 parsec survey (|Henrv et al.l[l99l . 

We use VB 10 as a standard to estimate the contrast 
ratio of the two components of GJ 164 in the R band 
where the astrometry data was taken. The 2MASS cat- 
alogue provides a (i? — K) = 6.8 color index for VB 
10, which combined with the dat a gives an apparent 
R = 16.6 magnitude for GJ 164B. iMonet et all (120031 ) 
give R = 12.4 for GJ 164A, which corresponds to a con- 
trast ratio c w 47. This gives an upper limit to the frac- 
tional light of the secondary: (3 — 0.021, which means 
that the hypothesis of a secondary of negligible visible 
luminosity would, at most, underestimate the mass ratio 
(cf. eq. H]) by 2.1 %. This systematic uncertainty needs 
to be added up to the statistical uncertainty deduced 
from the analysis performed in Section 13.21 

3. CHARACTERISTICS OF GJ 164 AB 

3.1. Orbital Parameters 

The aperture masking observations gathered in Table 
[T] provide the instant position of the companion to the 
primary. These observations alone suffice to constrain all 
seven orbital elements of the binary, independently of the 
parallax, with the semimajor axis expressed in angular 
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units. 

We use a m odel identical to t h e one described in Sec- 
tion 3.1.2 of iMartinache et al.l ()2007[ ) to fit a seven- 
parameter model to 14 observables. The observables are 
two coordinates for seven aperture masking data points. 
The seven orbital parameters are: the angular semimajor 
axis a, the eccentricity e, the longitude of the ascending 
node Wo, the inclination i, the argument of the periastron 
the period P and the epoch at periastron passage Tp. 
Our solution exhibits a final reduced xt = 0-61 for 7 de- 
grees of freedom. Note that the temporal coverage of 
these observations definitely rules o ut the possibi l ity of a 
~ 4 yr orbital period suggested bv iPravdo et al.l (|2004f ). 
The seven orbital parameters and the associated uncer- 
tainties are presented in Table ID 

We then proceed to a joined fit AO -I- astrometry, that 
is a thirteen-parameter fit to 40 observables. The extra 
parameters included in this new analysis are: the two 
components of the proper motion, the relative parallax 
n as well as the semimajor axis of the motion of the pri- 
mary around the center of mass of the system ai . With a 
2 mas uncertainty on t he individual astrom etric measure- 
ments, as assumed in iPravdo et al.l ()2004[ ). the solution 
to this global fit exhibits a final reduced xl = 0-73 for 
27 degrees of freedom. We therefore assumed astromet- 
ric errors of 1.4 mas so that the reduced xt is exactly 1. 
Table |4] compares the seven orbital parameters and their 
uncertainties obtained from the fit of the AO data only 
and the joined fit: both approaches produce compatible 
solutions. Figure [2] represents the orbit resulting from 
this joined fit, onto which are superposed the AO and 
astrometry data. 

3.2. Dynamical Masses 

Like lPravdo et al.l (|2004( l. we use a 2±1 mas correction 
to convert the relative parallax 11 = 70.0 ± 0.8 mas, to 
the absolute parallax: 11 — 72.0 ± 1.2 mas. We can 
now determine the actual relative semimajor axis of the 
binary: a = 1.1 ± 0.2 AU, as well as the total mass: 

MT = ayP^, (5) 

that is Mt = 0.343 ± 0.026 Mq. In the hmiting case 
where the luminosity of the secondary is negligible, the 
fractional mass of the secondary (cf. eq. [4]) is given 
by the ratio of the photocentric and relative orbit semi- 
major axis, that is f — 0.250 ± 0.010. It is remarkable 
that the fractional mass should be a ratio of two angular 
radii, for this makes this quantity parallax-independent. 
Section [2?2l shows there is at least a 3.8 magnitude differ- 
ence between the two components of the binary, which 
translates into a fractional light (3 — 0.021. The frac- 
tional mass / therefore requires non symetric uncertain- 
ties: / = 0.250j:°:S]3J. 

From the total mass and the fractional mass, one can 
infer the masses of both components: A/2 — 0.086]']q qq7 

and Ml — 0.257tQQ22- dynamical characteristics of 
the GJ 164 system are given in Table [6l 

4. SPECTROSCOPY 

Careful characterization of the atmosphere and deter- 
mination of abundances is necessary to understand the 
location of GJ 164 AB in a mass-luminosity (M/L) dia- 
gram, relative to an observational M/L relation or to nu- 



TABLE 4 
ORBITAL ELEMENTS 



Parameter 


AO only 


AO 4 


- STEPS 


a (mas) 


80.4 


± 


1.3 


80.5 


± 


1.2 


e 


0.157 


± 


0.012 


0.161 


± 


0.012 


i (deg) 


121.9 




0.9 


121.9 


± 


0.8 


fli (deg) 


272.4 


± 


1.3 


271.8 


± 


1.2 


LJo (deg) 


314.7 


± 


5.0 


311.7 




3.7 


P (days) 


734.3 




4.6 


736.9 


± 


1.7 


Tp (reduced JD) 


1868 


± 


19 


1856 
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TABLE 5 
ASTROMETRIC ELEMENTS 



Parameter Value 



Proper motion R.A. (mas) 


-324.5 


± 


0.3 


Proper motion Decl. (mas) 


-808.3 


± 


0.3 


Relative Parallax (mas) 


70.0 


± 


0.8 


Phot, semimaj. axis (mas) 


20.1 


± 


0.8 


Total Mass (M^) 


0.343 


± 


0.026 



TABLE 6 

DYNAMICAL MASSES: obtained from the 

JOINED ANALYSIS AO + STEPS. SEE TEXT FOR 
MORE DETAILS ON UNCERTAINTIES! 



Parameter 




Value 






Total Mass 


Mt 


= 0.343 ± 


0.026 


M@ 


Primary Mass 


Ml 


= 0.257 ± 


0.020 


Mq 


Secondary Mass 


M2 


= 0.086 ± 


0.007 


Me 


Mass ratio 




= 0.250 ± 


0.010 





merical simulations of M-Dwarf atmospheres. Determin- 
ing the overall metallicity of mid-to-late lype M-Dwarfs is 
however a delicate task. Indeed, with (V-K) > 5 for spec- 
tral types later than M-4, optical spectroscopy doesn't 
appear as an efficient approach, while in the infrared, 
there is no well-established proxy for metallicity that cal- 
ibrates with its optical counterpart. 

One strategy to solve this problem is to observe binary 
systems where the primary is a Solar-type star (F,G or 
K-Dwarf) of known metallicity and the secondary a M- 
Dwarf. Assuming that both components formed from 
the same original molecular cloud, they should share 
the same metallicity. In the context of addressing the 
planct-mctallicity correlatio n for M-Dwarf p lanet surveys 
lEdclstcin ct al. 2007; Muir head et al.ll2008l ) we have be- 
gun a separate study of a population of such objects 
(iRojas-Avala fc Llovdl l2008f ) . Initial results from this 
study and its conclusions for the metallicity of GJ 164 
are presented in this section. 

4.1. Observations and Data Reduction 

Our spectroscopic sample consists of five M4-5-Dwarfs 
(in complement to GJ 164 AB) associated to So- 
lar type stars, whose metallicity was measured by 
IValenti fc Fischeil ()2005[ ). Their designation, spectral 
type, V and K apparent magnitudes, as well as metallic- 
ities are listed in Table [T] 

Near-infrared spectra of these objects were obtained 
with the recently commissioned Trip leSpec spectrograph 
on the Palomar Hale Telescope (IWilson et al.l 120041 : 
iHerter et aLll2008f) . TripleSpec at Palomar is a A/ A A 
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1.20 



1.25 



1.30 



1.35 



Fig. 3. — J-band spectra of the stars ordered from top to bottom by decreasing metallicity (listed in Table (7)1. The most prominent 
spectral features are highlighted. 




X (|im) 



Fig. 4. — K-band spectra of the stars ordered from top to bottom by decreasing metallicity (listed in Table The most prominent 
spectral features are highlited. 
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Fig. 5. — Equivalent widths of some of the most prominent 
lines in the sample spectra with 1 a errors versus (V-K) colors. A 
diamond represents GJ 164 AB, metal-rich objects ([M/H] > 0.0) 
are represented by asterisks and metal-poor objects ([M/H] < 0.0) 
by crosses. 
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Fig. 6. — Spectrum of GJ 164 (top) compared to the PHOENIX 
models oflHauschildt et al.l I I1999I) for Solar metallicity: the models 
(labeled PHOENIX: M4) reproduce most features of the spectrum. 
The presence of a M8-Dwarf companion with a 5:1 K-band contrast 
ratio (labeled PHOENIX:M4-f M8) only modifies the strength of the 
GO bands, as shown in the plot of the difference between the two 
models (bottom). 



TABLE 7 
The Sample 



Name 


Sp. Typ. 


V 


K 


[M/H] 






[mag] 


[mag] 


dex 


GJ 324 B 


M4 


13.15 


7.67 


4-0.25 ± 0.03 


GJ 544 B 


M4 


14.50 


9.59 


-0.15 ± 0.03 


GJ 611 B 


M4 


14.20 


9.16 


-0.49 ± 0.03 


GJ 166 G 


M4.5 


11.17 


5.96 


-0.08 ± 0.03 


NLTT 25869 


M5 


14.50 


8.64 


-1-0.27 ± 0.03 



« 2500-2700 cross-dispersed near-infrared spectrograph 
with a broad wavelength coverage across 5 simultaneous 
orders (1.0-2.4 fim) in echelle format. Its entrance slit is 
Ix 30 arcseconds and a notable feature is that it has no 
moving parts. 

The data were reduced with an IDL-based data reduc- 
tion pipeline developed by P. JVIuirhead for TripleSpec 
at Palomar. The data were sky-subtracted using a sky- 
frame made by median combination of the 5 exposures 
on different positions along the slit of the object. To 
correct for telluric absorption features and flux-calibrate 
their spectra, an AO V star was observed as close to 
the science object star airmass as possible. Each sky- 
subtracted exposure was then divided by a normalized 
flat-field, wavelength calibrated and fully extracted. The 
spectra are flux-calibrated and telluric corrected using 
the IDL- based code x t ellco r_general described in the 
paper bv 'Vacca et all (|2003l ). 

The final J- and K-band spectra along with the most 
prominent features are shown in Figs. [3]and[4l respec- 
tively. 

4.2. Equivalent Widths 

The equivalent width (EW) of six prominent features 
in the J and K-band spectra were estimated using the the 
IDL-based function measure_ew by N. Konidaris and J. 
Harker. Each spectrum was normalized and the equiv- 
alent widths were calculated by computing the ratio of 
the area of a feature to a pseudocontinuum. This pseu- 
docontinuum is a linear interpolation of clean regions on 
either side of the feature. Uncertaities in the equivalent 
width were obtained usin g the procedure described in 
iSembach fc Savagd ([199^ . 

The equivalent widths of K I (1.178 fim), Fe I (1.189 
^m) and Mn I (1.290 ^m) in J-band as well as Na I (2.206 
/zm & 2.209 ^m), Ca I (2.261 /im, 2.263 /im & 2.265 nm) 
and 12CO(2,0) (2.293 ^m) in K-band are shown in Fig. 
[51 Crosses depict the stars with [M/H] < [M/H]© (GJ 
544 B, GJ 611 B and GJ 166 C) and asterisks depict 
the stars with [M/H] > [M/H]© (GJ 324 B and NLTT 
25869). A diamond represents GJ 164 AB. 

4.3. Metallicity of GJ 164 

Except for the Fe I lines tha t do not exhib it any 
correlation with spectral type (jCushing et al.l [2005), 

most of the spectral features show n in Fig . [5] 

are strongly tem perature-dependent ()Ali et al.l 119951 : 
iJones et alJil996l ). The dispersion in EW between ob- 
jects with roughly the same temperature and gravity will 
therefore provide insights on the metal abundances, and 
the spectrum of GJ 164 shows features of strength com- 
parable to the metal-rich stars of our sample, especially 
in K-band (cf. Fig. [4]). Except for K I (1.178 ^ni) and 
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Fig. 7. — GJ 164 A and B in a Mass-Luminosity diagram. The 
absolute K magnitudes are compared to the BCAH98 isochrones 
for 0.1 and 1.0 Gyr. Asymmetric confidence intervals on the masses 
of both components are due to the uncertainty on the contrast ratio 
at the wavelength where the astrometry was taken. 

Fe I (1.189 //m), the EW of its lines (cf. Fig. [5]) make 
GJ 164 more likely to be metal-rich than metal-poor. 

The high value of GJ 164's ^^C0{2fl) (2.293 fim) EW 
(Fig. [5l bottom right panel) can be attributed to its sec- 
ondary. Indeed, the data presented in Section 12.11 shows 
that the contrast ratio between GJ 164 A (M4) and GJ 
164B (M8) is close to 5:1 in K. Since the strength of 
the CO bands increases with later types in the M-Dwarf 
sequence, GJ 164 B (unresolved with TripleSpec) adds 
a non-negligible contribution to the strength of these 
molecular features. Figure [5] compares the K-band spec- 
truni of GJ 164 to PHOENIX spectra (jHauschildt et all 
119991 ). A M4-Dwarf spectra alone as well as a spectrum 
composed of a M4 and a MB with a 5:1 contrast ratio, 
aU for [M/H]=0.0. Although the strength of the Na I 
doublet and the Ca I triplet in the PHOENIX spectra dif- 
fer from our data, the strengths and shapes of the CO 
bands predicted by the models match the spectra of GJ 
164. The difference between these two spectra at the 
bottom of Fig. [H] shows that the only features affected 
by the MS companion at these wavelengths arc the CO 
bands. 

5. MASS-LUMINOSITY RELATIONS 

The individual apparent magnitudes calculated in Sec- 
tion 12. 2i combined with the absolute parallax deduced 
from the dynamical analysis performed in Section 13.21 
determine the following absolute H and K magnitudes: 
Mh = 7.721 ±0.051 and Mk = 7.407 ±0.047 for GJ 164 
A, Mh = 9.557 ± 0.051 and Mk = 9.128 ± 0.047 for GJ 
164 B. 

The (H-K) colors, corroborated by the infrared spec- 
troscopy data presented in Section [3] indicate that GJ 
164 is at least of Solar metallicity. Figure [7] compares 
the location of GJ 164 A and B in a Mass- luminosity di- 
agram, relative to the low-mass Solar metallicity models 



of lBaraffe etldl (|1998D (hereafter refered to as BCAH98) 
for the K-band. For masses below 0.4 Mq (but above the 
substellar limit), models predict very little evolution be- 
tween 0.5 and 5 Gyr. Fig. [7] therefore only plots two 
isochrones, for 0.1 and 1 Gyr. 

The aperture masking data places the primary less 
than l-tj away from the 1 Gyr isochrone, and 3-tT away 
in mass from the 0.1 Gyr isochrone. The hypothesis of 
an evolved age is supported by the kinematics of the sys- 
tem. Indeed, with the proper motion (-324.5, -808.4) mas 
yr^^, the absolute parallax 11 = 72 ± 1.2 mas derived in 
Section [321 as w ell as the radi a l velo city V ~ —29.9 km 
s^^ measured bv lPravdo et al.l (|2004f ). one can calculate 
the Galactic space velocity {U, V, W) = (-24, -20, -47) 
km after correction for standard solar motion. 
When compared to the properties of the nearest young 
movin g groups (cf. Table 1 of iLopez-Santiago et al.l 
(|2006f )). GJ 164 exhibits a velocity component perpen- 
dicular to the Galactic plane too large to be associated 
to any of these groups. Fro m the kinemat ic properties of 
the M-Dwarf survey of Rei d et al.l ()2002D . it is however 
impossible to conclude whether GJ 164 belongs to the 
thin or the thick disk. 

The case of GJ 164 B is a little bit more puzzling. In- 
deed, as the mass decreases, the luminosity predicted by 
the models becomes a very steep function of the mass. 
GJ 164 B lies 1.5 a away from the 1 Gyr isochrone hor- 
izontally, which is only marginally in conflict with the 
models. However, based on these dynamical mass and 
luminosity measurements, the models underestimate the 
luminosity of a <0.1 Mq star by roughly one magnitude. 
A way to accomodate for this discrepancy would be for 
the GJ 164 A to be a tight binary, like for exam ple the 
case of GJ 802 exposed i n [Ireland et al.l (l2008f) . Spec- 
troscopic observations by iPravdo et al.l ()2004f ) however 
showed no evidence that GJ 164 A is a spectroscopic 
binary. 

To provide a stronger case. Figure [5] plots more of the 
BCAH98 isochrones (0.1, 0.2, 0.3, 0.5 and 1 Gyr), this 
time in the plane mass ratio f (cf. eq. H]), against the K- 
band magnitude difference. AK is a direct measurement 
made with the aperture masking intcrferometry data (cf. 
Table [2]) while f is a derived product of the joined fit, 
obtained in Section 13.21 Both observables are parallax- 
independent, which makes them more robust than the 
individual masses and luminosities. 

Although taking into account all the uncertainties of 
our analysis by generating families of models within ± 1 
standard deviation of the total mass of the binary, the 
1.5-cr discrepancy between the data and the models pre- 
sented in Fig. [7] now turns out to be over 5 a away from 
the 1 Gyr isochrone: the models presented in Fig. [5] 
clearly favor a very young age (100 Myr), which accomo- 
dates for the apparent excessive luminosity of GJ 164 B, 
but is incompatible with the kinematics of the system. 
Note that the same conclusion holds when performing a 
similar analysis using H-band mass-luminosity relations. 

6. CONCLUSION 

^ The sign convention is the one of the IDL astrolib gaLuvw pro- 
cedure, with U positive toward the Galactic anti-center, V positive 
in the direction of Galactic rotation, and W positive toward the 
North Galactic Pole. 
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Fic;. 8. The mass ratio and K-band contrast ratio of GJ 164 
(datapoint of coordinates (0.25, 1.8)) arc compared to BCAH98 
isochrones for ages between 0.1 and 1 Gyr. Each age generates a 
family of models taking into account the uncertainty on the total 
mass of the binary. 

GJ 164 is an astromctric binary, whose orbit was re- 
solved by aperture masking interferometry observations 
at Palornar and Keck. The consistency of the data is 
excellent, and exhibits a 1-2 mas precision in average for 
each measurement, made below what is usually accepted 
as the resolution limit of a telescope. This data, com- 
bined with earlier astrometric measurements, provides 
the following dynamical masses: 0.257 ± 0.020 M0 for 
the primary and 0.086 ± 0.007 Mq for the secondary. 

Analysis of its infrared spectrum reveals that the M- 
Dwarf binary GJ 164 is at least of Solar metallicity. 
When attempting to use theoretical models for Solar 
metallicity very low-mass stars to derive an age from 
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the mass and luminosity of both components, we find 
that the models do not adequately fit the data, requir- 
ing very young ages to accomodate for an overluminous 
secondary. The models therefore underestimate the lu- 
minosity of very low-mass stars that have settled on the 
main sequence. 

The precision of the dynamical data we present here, 
of the order of 10 percents, however does not provide a 
strong constraint for the models. A precise measurement 
of the contrast ratio in the R band, where the astrome- 
try data was taken would already help to better constrain 
the mass ratio and ultimately improve the constraint on 
the individual masses, and make our statement about the 
models stronger. The very small separation (at most, 90 
mas) makes this measurement difficult, even for HST. 
But because it lies so close from the substellar limit, GJ 
164 B is a target that deserves most attention. In the 
meantime, a Radial Velocity curve would provide an in- 
dependent measurement of the mass ratio. With an ex- 
pected Radial Velocity amplitude of 2.8 km/sec, and a 
period of two years, observations of GJ 164 with Triple- 
Spec and TEDI will settle this uncertainty. 
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